Isospin effects on sub-threshold kaon production at intermediate energies 
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We show that in collisions with neutron rich heavy ions at energies around the production thresh- 
old K° and K + yields might probe the isospin dependent part of the nuclear Equation of State 
(EoS) at high baryon densities. In particular we sugg est the K°/K+ ratio as a promising observ- 
able. Results obtained in a fully covariant relativistic transport approach are presented for central 
Au + Au collisions in the beam energy range 0.8— 1.8 AGeV. The focus is put on the EoS influence 
which goes beyond the collision — cascade picture. The isovector part of the in-medium interac- 
tion affects the kaon multiplicities via two mechanisms: i) a symmetry potential effect, i.e. a larger 
neutron repulsion in n-rich systems (isospin fractionation); ii) a threshold effect, due to the change 
in the self-energies of the particles involved in inelastic processes. Genuine relativistic contributions 
are revealed, that could allow to directly "measure" the Lorentz structure of the effective isovector 
interaction. 
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For at least two decades particle production has been 
suggested as a useful tool to constrain the poorly known 
high density behaviour of the nuclear equation of state 
(EoS) P, Q • This information is crucial both to learn 
about the effective nuclear interaction and to get a 
deeper insight into many interesting astrophysical phe- 
nomena. In particular pion and (subthreshold) kaon pro- 
ductions have been extensively investigated both theoret- 
ically 0, H H H H 0, 11 and experimentally @, G3 El, 
leading to the conclusion of a soft behaviour of the 
EoS at high densities, see the recent refs. ^2 liflj - 
Kaons (K°,K + ) appear as particularly sensitive probes 
since they are produced in the high density phase al- 
most without subsequent reabsorption effects. At vari- 
ance, antikaons (K°,K~) are strongly coupled to the 
hadronic medium through strangeness exchange reac- 
tions 0, Q, E| . In this paper we show that the isospin 
dependence of the K°' + production cross sections can be 
also used to probe the isovector part of the EoS: we 
propose the K°/K + yield ratio as a good observable to 
constrain the hi gh density behavior of the symmetry en- 
ergy, E sym , dill. 

In heavy ion collisions (HICs) at intermediate ener- 
gies (1-2 AGeV) a transient state of highly compressed 
and heated matter is created, with density up to 3-4 
times the saturation value. Several dynamical observ- 
ables have hence been proposed to characterize the den- 
sity dependence of E sym , such as charge ratios and col- 
lective isospin flows for highly energetic nucleons and pi- 
ons 0, 0, EH El ES El and isospin transparency as 
expressed in transport ratios 0j EH El< Controver- 
sial deductions arise from different relativistic and non- 
relativistic models, see ^3- 



We show that, within a covariant description of nu- 
clear dynamics, from the K° / K + ratios we can directly 
investigate the Lorentz structure, i.e. the scalar-vector 
decomposition, of the isovector sector of the effective 
in-medium hadron Lagrangian. Some promising indica- 
tions have been recently obtained in nuclear matter cal- 
culations E3> here we present results for realistic open 
systems, i.e. for collisions of neutron-rich heavy ions in 
the energy range around the kaon production threshold 
(1.56 AGeV). 

Using a transport model, derived within the Rela- 
tivistic Mean Fiel d ap proximation (RMF) of Quantum- 
Hadro-Dynamics [28j, we analyze pion and kaon pro- 
duction in central 197 Au + 197 Au collisions in the 0.8 — 
1.8 AGeV beam energy range, with different RMF ef- 



fective field choices for E* 



We will compare results 



of three Lagrangians with constant nucleon-meson cou- 
plings (NL... type, see 0, Eof) and one with density 
dependent couplings (DDF, see |2fj), recently suggested 
for apparent better nucleonic properties of neutron stars 
29] . In order to isolate the sensitivity to the isovector 
components we use models showing the same "soft" EoS 
for symmetric matter El . This is achieved via Non- 
Linear (NL...) contributions of the isoscalar scalar a and 
vector lo mesons in the constant coupling cases, and via 
suitable density dependences in the DDF model. 

For the isovector part, in the simple NL choice the 
mean field is not isospin dependent and the symmetry 
energy is just due to kinetic (Fermi) contributions. The 
NLp model contains an isovector-vector effective field 
(p meson), which leads to a splitting of the vector self- 
energies between protons and neutrons, and to a more 
repulsive force experienced by neutrons with respect to 



protons in neutron-rich matter [l?J. In the NLpS model 
also an isovector-scalar S field is included, which gives 
an effective mass splitting between protons and neutrons 

nam. 

With both p- and (5-meson couplings, f Pj g, E sym can 
be written as [3l| 
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leading to a partial cancellation of scalar and vector 
components, equivalent to that of the er and u> fields in 
the isoscalar sector. Since the 5 field couples to the scalar 
density, f$ is multiplied with a density dependent quench- 
ing factor in Eq.£[J. It is then clear that the increase by 
nearly a factor of three in f p in the NLpS model is nec- 
essary to compensate the attraction due to the scalar S 
field, in order to get the same bulk asymmetry parame- 
ter a 4 = 30.5AjfeV^ at saturation, [2fJ,|3j|. Consequently 
the inclusion of a 5 field leads to a stiffer symmetry en- 
ergy at high baryon densities |s3, an d to larger vector 
self-energies for nucleons, as discussed in |lJl|22j]. In the 
DDF model the f p is exponentially decreasing with den- 
sity, resulting in a rather "soft" symmetry term at high 
density [HEi^. 

We like to note the genuine relativistic nature of these 
mechanisms that will directly influence transport dynam- 
ics and particle production. A transport code with only 
relativistic kinematics combined to density dependent in- 
teractions will miss these effects. 

In the energy range considered here the nucleon- 
nucleon inelastic channels can be restricted to the exci- 
tation of the lowest mass resonance A(1232) and pertur- 
bative kaon (K +, °) production through baryon-baryon 
collisions BB — ► BYK, where B stands for nucleons or 
resonances and Y for hyperons (A, I]^ ). Pions are pro- 
duced via the decay of the A(1232) resonance and - after 
propagation and rescattering - can contribute to the kaon 
yield through collisions with baryons: ttB — ► YK. All 
these processes are treated within a relativistic hadronic 
transport model of Boltzmann-Uehling-Uhlenbeck type 
(RBUU) ] i. e. including a hadron mean field propaga- 
tion The latter point, which goes beyond the 
"collision cascade" picture, is essential for particle pro- 
duction yields since it directly affects the energy balance 
of the inelastic channels, as extensively analysed in the 
following. Further details can be found in ref. , where 
all the parametrizations used for the corresponding cross 
sections are also specified. 

Due to the K +fi long mean free path their in-medium 
widths can be expected to be small and an on-shcll quasi- 
particle treatment appears suitable. Here, as in [2j|, we 
do not generally use potentials for kaons and propagate 
them as free particles. This simplification is justifiable, 
since our main aim here is to show how the Lorentz struc- 
ture of the symmetry energy influences the K° / K + ratio, 
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FIG. 1: Time evolution of the A ±,0 ' ++ resonances (left) and 
the pions 71-*'° (right) for a central (6 = fm impact param- 
eter) Au+Au collision at 1 AGeV incident energy. Transport 
calculation using the NL, NLp, NLpS and DDF models for 
the iso-vector part of the nuclear EoS are shown. 



which should not sensitively depend on the kaon-nucleon 
potential as this appears roughly the same (slightly re- 
pulsive) for both K° and K+ 

Fig. ^ reports the temporal evolution of A ±,0,++ res- 
onances and pions (n^' ) (left panel) and kaons (AT +, °) 
(right panel) for central Au+Au collisions at 1 AGeV. 
Results are shown for all the tested models. 

It is clear that, while the pion yield freezes out at times 
of the order of 50 fm/c, i.e. at the final stage of the 
reaction (and at low densities), kaon production occurs 
within the very early stage of the reaction, and the yield 
saturates at around 20/m/c. Kaons are then suitable to 
probe the high density phase of nuclear matter. This is 
not the case for pions, which suffer from reabsorption and 
isospin exchange processes that modify both the absolute 
primordial yield and the ir~ /ir + ratio. In fact from Fig. 
^we see that the results for different models are rather 
similar, and thus pion multiplicities depend only weakly 
on the isospin part of the nuclear mean field. However, 
a slight increase (decrease) in the 7r~ (7r + ) multiplicity is 
observed when going from the NL to the NLp and then 
to the NLpS model, i.e. increasing the vector contribu- 
tion f p in the isovector channel. This trend is more pro- 
nounced for kaons, see Fig. fright panel, due to the high 
density selection of the source and the proximity to the 
production threshold. The results for the DDF model, 
density dependent couplings with a large f p decrease at 
high density, are fully consistent. They are always closer 
to the NL case (without isovector interactions) but the 
difference still seen for is completely disappearing 

for K°' + , selectively produced at high densities. 

When isovector fields are included the symmetry po- 
tential energy in neutron-rich matter is repulsive for neu- 
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trons and attractive for protons. In a HIC this leads to a 
fast, pre-equilibrium, emission of neutrons towards more 
symmetric nuclear matter. Such a mean field mecha- 
nism, which has been extensively investigated, often re- 
ferred to as isospin fractionation 0, , is responsible 
for a reduction of the neutron to proton ratio during the 
high density phase, with direct consequences on particle 
production in inelastic NN collisions. 

Threshold effects represent a more subtle question. 
The energy conservation in a hadron collision in general 
has to be formulated in terms of the canonical momenta, 
i.e. for a reaction 1 + 2 — > 3 + 4 as 

s in = (&i + &2 ) 2 = (^3 ^4 ) 2 = s out (2) 

Since hadrons are propagating with effective (kinetic) 
momenta and masses, an equivalent equation should be 
formulated starting from the effective in-medium quanti- 
ties k*^ = k^ - T,^ and m* = m + £ s , where S s and W 
are the scalar and vector self-energies, respectively, de- 
pending on the isovector channel structure. In particular, 
for the general aujpS case one obtains for the self-energies 
of protons and neutrons: 

£»(p,n) = -foPs ± fsPs3 (3) 
E"(p,n) = UFTfpjg, ■ (4) 

(upper signs for neutrons), where p s — p sp + p sn , j a = 

3p + 3n,Ps3 = Psp ~ Psn, j% = 3p ~ 3n are tlie total and 

isospin scalar densities and currents and f a ,u,p,s are the 
coupling constants of the various mesonic fields, |3^. 

In reactions where nucleon resonances, especially the 
different isospin states of the A resonance, and hyperons 
enter, also their self energies are relevant for energy con- 
servation. We specify them in the usual way according 
to the light quark content and with appropriate Clebsch- 
Gordon coefficients |27j . 

In the most general case the isovector scalar and vector 
self-energies enter the new threshold condition for a given 
inelastic process 

s in > (m* 3 + £° + m\ + £°) 2 - (E 3 + S 4 ) 2 . (5) 

The condition of energy conservation in inelastic 
hadron collisions will influence the particle production 
in two different ways. On one hand it will directly deter- 
mine the thresholds and thus the multiplicities of a cer- 
tain type of particles, in particular of the sub-threshold 
ones, as here for the kaons. Secondly it may favour or 
penalize reactions, because the self-energies in the final 
channel are more attractive or repulsive than in the ini- 
tial one, and consequently the phase space in the final 
channel is larger or smaller. 

Hence, we clearly see that the competition between 
scalar and vector isovector fields is responsible for the 
isospin fractionation and modifies particle production 
rates. In particular while the scalar and vector isovector 



fields tend to cancel in the symmetry term, see Eq.lQJ, 
they can have very different dynamical effects, as already 
noted in the flow analysis of ref . [lflj . In fact mean field 
and threshold effects are acting in opposite directions on 
particle production and might compensate each other. 

As an example, nn collisions excite A~'° resonances 
which decay mainly to n~ . In a neutron-rich matter the 
mean field effect pushes out neutrons making the matter 
more symmetric and thus decreasing the ir~ yield. The 
threshold effect on the other hand is increasing the rate 
of 7r -, s due to the enhanced production of the A - res- 
onances: now the nn — > p/Sr process is favored (with 
respect to pp — > nA ++ ) since more effectively a neutron 
is converted into a proton. 

We have to note that such interplay between the two 
mechanisms cannot be fully included in a non-relativistic 
dynamics, in particular in calculations where the baryon 
symmetry potential is treated classically |2l|, |22, |23j . 
A typical example is the strength of the isovector-vector 
p-couphng which is linked to the symmetry energy but it 
is largely varying with the Lorentz structure of the isovec- 
tor interaction. Actually aim of this paper is to show 
that the data could be used to directly probe the covari- 
ant structure of the effective interaction in the isovector 
channel at high baryon densities. 

In the left panel of Fig. \l\we see that A - (and tt~) pro- 
duction increases when going in the direction of larger f p 
couplings NL — > DDF — > NLp — ■> NLpd. As discussed 
before, the isospin fractionation lowers the n/p ratio and 
thus decreases the A~ (and 7r~) yield, while the threshold 
mechanism increases it due to the larger nn self-energies. 
In fact the latter effect appears to overcompensate the 
former. The opposite is true for the A +,++ (and 7r + ) 
production. 

As seen in the right panel of Fig. ^ a similar argument 
holds for K° and K + mesons, which mainly come from 
nn (or ir~n) and pp (or 7r + p) collisions, respectively, and 
thus exhibit the same trend as tt~ and tt + . However, the 
isospin effect in this case is more pronounced because the 
changes in the self energies for the different models play a 
more crucial role close to the kaon production threshold. 

Finally the beam energy dependence of the tt^/tt + 
(left) and K° / K + (right) ratios is shown in Fig. [3 At 
each energy we see an increase of the yield ratios with 
the models NL -> DDF -> NLp -> NLpS. The effect 
is larger for the K°/K + compared to the tt~ /tt + ratio. 
This is due to the subthreshold production and to the 
fact that the isospin effect enters twice in the two-steps 
production of kaons, see j3j|. Between the two extreme 
DDF and NLpS isovector interaction models, the varia- 
tions in the ratios are of the order of 14 — 16% for kaons, 
while they reduce to about 8 — 10% for pions. Interest- 
ingly the Iso- EoS effect for pions is increasing at lower 
energies, when approaching the production threshold. 

We have to note that in a previous study of kaon pro- 
duction in excited nuclear matter the dependence of the 
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K°/K + yield ratio on the effective isovector interaction 
appears much larger, about ten times more for a system 
with the Au asymmetry (see Fig. 8 of ref.(27|). The point 
is that in the non-equilibrium case of a heavy ion collision 
the asymmetry of the source where kaons are produced 
is in fact reduced by the n — > p "transformation" , due to 
the favored nn — » pA~ processes. This effect is almost 
absent at equilibrium due to the inverse transitions, see 
Fig. (3) of ref.[i3|. Moreover in infinite nuclear matter 
even the fast neutron emission is not present. 

Thus in an open system the signal appears weaker but 
still accessible experimentally, if it can be shown to be 
robust. In order to further stress the distinction between 
effects of the stiffness of the symmetry energy and the 
detailed Lorentz structure of the isovector part of the 
effective Lagrangian, we also show the results for the 
K°/K + with another parametrization of E sym . This 
model, NLDDp, is a variant of NLp with a density 
dependent p-coupling, built in such a way as to repro- 
duce the same stiffer E sym {pB) of the NLpS model (see 
also ref. 0). The results for the n~ /tt+ and K°/K+ 
ratios are shown in Fig. [21 for Eb eam — 1.0 AGeV as 
triangles. We see that they are closer to the NLp re- 
sults (with a constant f p ) than to the ones of the NLpS 
choice which has the same iso-stiffness. This nicely con- 
firms that the differences observed going from the NLp 
to the NLpS parametrization are not due to the slightly 
increased stiffness of E sym (pB), but more specifically to 
the competition between the attractive scalar <5-field and 
the repulsive vector p-field in the isovector channel, which 
leads to the increase of the vector coupling (see the com- 
ments to the Eq.JD). Hence the differences between the 
vector self-energies of neutrons and A - resonances vs. 
the proton and A ++ ones are greater in the NLpS model 
than in other models without scalar isovector fields and 
the production cross sections are consequently affected. 
As a result, the subthreshold K° (K + ) yield is sensitively 
enhanced (reduced) in the NLpS case. 

In the same Fig. [21 we also report results at 1.2 AGeV 
for the 132 Sn + 124 Sn reaction, induced by a radioac- 
tive beam, with an overall larger asymmetry (open sym- 
bols). The isospin effects are clearly enhanced. Finally 
we have performed some calculations including kaon po- 
tentials, derived from a chiral lagrangian as suggested in 
The K° /K + ratio appears not very sensitive to this 
in-medium modification of the kaon single particle ener- 
gies. This is somewhat expected since we have a similar 
small repulsive contribution for both A'° :+ mesons. 

From the present discussion we conclude that sub- 
threshold kaon production could provide a promising tool 
to extract information on the isovector part of the nuclear 
interaction at high baryon density. We have seen that, at 
beam energies below and around the kinematical thresh- 
old, the K°/K + inclusive yield ratio is more sensitive to 
the Lorentz structure of E sym than the 7r~/7r + . 

We would like to stress the two most important re- 
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FIG. 2: iv~ /tt + (upper) and K + /K° (lower) ratios as a func- 
tion of the incident energy for the same reaction and models 
as in Fig. In addition we present, for Et ea m = 1 AGeV, 
NLp results with a density dependent p-coupling (triangles), 
see text. The open symbols at 1.2 AGeV show the corre- 
sponding results for a 132 Sn -f 124 Sn collision, more neutron 
rich. Note the different scale for the ir~ /-k + ratios. 



suits of our study: i) We have shown that isospin effects 
are important not only at the mean-field level (isospin 
fractionation) but they also influence significantly parti- 
cle production cross sections. As a matter of fact, we 
observe that, for the reactions studied here, the modi- 
fications induced in the inelastic vertices represent the 
dominant effects, ii) At relativistic energies, due to the 
Lorentz structure of the isovector nuclear interaction, the 
isotopic content of particle emission is not directly related 
to the symmetry energy value, but it can be rather con- 
sidered as a measure of the strength of isovector vector 
channel. 



We note that the isospin effects on the kaon inclusive 
yield ratios at the freeze-out appear not too strong, al- 
though experimentally accessible. It seems important to 
select more exclusive kaon observables, in particular with 
a trigger related to an early time A-production. A trans- 
verse momentum selection of pion yields, corresponding 
to a higher density source, should also be rather sensitive 
to isospin effects, in particular at lower energies, closer to 
the production threshold. A large asymmetry of the col- 
liding matter is in any case of relevance. In this sense our 
work strongly supports the study of particle production 
at the new relativistic radioactive beam facilities. 
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